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1 Introduction 

Geometric phases have received a great deal of attention since their descrip- 
tion by Berry Jl| . The reasons are clear. They are a fundamental property of 
many quantum mechanical systems. They also have a beautiful description 
in terms of differential geometry and fiber bundles || which is directly related 
to gauge field theory (see for example ||). Their physical importance was 
known long before the excitement about them in the mid 80's 0], ||, ||. In 
spite of all the attention, there have been few worked out examples and the 
examples that have been worked out, the descriptions haven't been straight- 
forward. The most well-known example is that of a two state system, namely 
a magnetic dipole, in a magnetic field. This was the original example given 
by Berry |TJ. Wilczek and Zee originally pointed out that there could exist 
non- abelian geometric phases JF|. Later people studied fermionic systems 
with a quadrapole Hamiltonian |J. 

Uhlmann later developed machinery, namely a parallel transport [HJ , for 



describing the non-abelian geometric phases associated with density matri- 
ces. However, this was never applied to three state systems. Arvind et al 
and Khanna et al studied the geometric phases for three state systems that 
involve pure state density matrices |TT|, |T2| with a parameterization that 
was somewhat ad hoc. Mostafazadeh looked at a way of calculating the non- 
abelian geometric phases for a three state system with a two-fold degeneracy 



[ 13fl also with ad hoc coordinates. These topics will be brought together here. 
In this paper the objective is to use explicit SU(3) representations to 
extend and/or simplify several aspects of three state systems. 

1. The expression for the density matrix for three state systems. 

2. The identification of the parameter spaces of these systems. 

3. The calculation of the abelian geometric phases for three state systems. 

4. The calculation of the non-abelian geometric phases of three state sys- 
tems with a two-fold degeneracy. 

An obvious example of a three state system would be a spin one particle 
in a magnetic field. If this external magnetic field is "slowly" rotating then 
we may have the conditions for an adiabatic change in phase. For a proper 
description of what "slowly" means in this context see jn[] . Unless otherwise 
stated, this paper will be concerned with the adiabatic geometric phases 
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Table 1: The Gell-Mann Matrices 

although with some work this could be extended to non-adiabatic phase 
changes. 

2 The Density Matrix for a System with Three 
Quantum States 

As is demonstrated in the next section the density matrix can be parameter- 
ized by the action of an SU(3) transformation. This will prove convenient 
for many calculations and is immediately generalizable to a system with an 
arbitrary number of states. See |TS|] and below for a discussion. 

The density matrix for general pure state three-level systems is given in 
[n| and | IHfl . It can be represented in the following way: Let ip be a state in 



a three dimensional complex Hilbert space Ti^. The density matrix is the 
matrix p described by (in analogy with two state systems): 

p = ^ = 1^^1 = 1(1 + ^3^ .X) (1) 

Here the dagger denotes the hermitian conjugate, n is a real eight dimensional 
unit vector, A represents the eight Gell-Mann matrices. The dot product is 
the ordinary sum over repeated indices n r X r . The (•, •) is the inner product 
on the space 7i^ . The pure state density matrix satisfies: 

pt = p 2 = p > Trp = L 



This is equivalent to the following conditions on n: 

n* = n n • n = 1 n-kn = n. (2) 

The star product is defined by 

(a -k b)i = V3d ijk ajb k (3) 

where the dijk are the components of the completely symmetric tensor ap- 
pearing in the anticommutation relations 

4 
{Aj, Xj\ = —dij + 2dijkXk- 

Explicitly the nonzero d^k are 

^118 — ^228 = ^338 = ~^888 = ~ 7= ^448 = ^558 = ^668 = ^778 = — r- 
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-doA7 — do^F, — d-tdd — d. 



146 — "157 — ~~ "247 — "256 — "344 — "355 — — "366 — — "377 



da 



3 Parameter Spaces for Three State Systems 

The parameter space of states of the three state systems can easily be seen 
to be coset spaces of SU(3). The Euler angle parameters are a particularly 
convenient way in which to see this. 

A representation of the coset space SU(3)/U(2) and of the density matrix 
for the pure states of a three state system may be obtained in terms of the 
Euler parameters given in fl6 |. There the group SU(3) is parameterized by 



D(a,(3, 7, 9, a, b, c, 



e (iXsa) e (iA2i9) e (iAs7) e (*A6fl) e (iX3o) e (iA 3 6) e (iX3c) e (iA8*)_ 



With this parameterization and the explicit representation of the correspond- 
ing adjoint representation in terms of the Euler angle parameters in [fL7| , a 



m 



parameterization of the density matrix of the three state system may be ob- 
tained by the following projection which is analogous to the Hopf map given 

-(1 - V3X 8 



x 



n(D) =D 



D-\ 



(4) 



Here x G SU(3)/U(2) and D represents a point in the space SU(3). This 
projection is clearly invariant under the right action of a U(2) operation 
defined by U E U(2) with 



U = e (i\ S a') e (i\2b') e (i\3c') e (i\ 8 <f>')_ 



This then defines the projection from 577(3) to SU(3)/U(2). Since the second 
term in equation (4) is simply an adjoint action on As, it can be read directly 

TT| . There the matrix Ra that satisfies 



from the equations given in 



l '.y 



DXiD' 



KijAj 



(5) 



was given explicitly. 

One may of course note that the projection operator is not unique. Any 
3x3 matrix with a one on its diagonal would be invariant under a U(2) 
subgroup and would represent a pure state. It is, however, rather convenient 
in this parameterization to use this particular matrix: 



1 



1- V3\ & 




so that it is clear that the upper left 2x2 matrix of zeros will be unaffected 
by an SU(2) transformation in that block. This matrix could be substituted 
for another that has one 1 on a diagonal and zeros elsewhere and still be 
invariant under (another) SU(2). The invariance of this with respect to an 
overall phase gives the U(2) invariance. 
Now equation (II) can be rewritten as 



x 



;i - VzRsjXj) 



-{l+VZnAj) 



(6) 



where we identify the R$j as the components of a vector that satisfies those 
properties given in equation ([!]). This can be viewed as an arbitrary rotation 
of the vector As with an adjoint action of the group (equation (|j)) and of 
course it is now clear that x is identified with p. 
The vector n has the following components. 



n-i 



R 



M 



V3 
cos 2a sin 2/5 sin 2 9 

2 H 



n 2 = — i?82 = — sin 2a sin 2(3 sin 2 9 
n 3 = —Rs3 = — cos 2(3 sin 2 9 

n 4 = — i? 84 = -— cos(a + 7) cos/3sin26 l 

/3 
n 5 = — i? 85 = — — sin(a + 7) cos/3sin26 l 

\/3 

n 6 = — i? 86 = cos (a — 7) sin /3 sin 2$ 

\/3 

727 = — R 8 j = — — sin (a — 7) sin (3 sin 29 

3 

n s = -R 88 = -l + -sm 2 9 (7) 

From this, using the equations n, = ip^iip, h follows that 

e i(a+7) cos /3 sin^ 
ijj = e ix \ e-^-T) sin /5 sin^ |. (8) 

cos^ 

This may be recognized as the third column of the S'L r (3) matrix I? above, 



thus agreeing with the calculation given in [1 Ij and coming full circle in the 
analysis. In this case the overall phase \ ma y be identified as — 20/v^3 in 
the matrix D. In section [| it will become clear why this works and it will be 
generalized for the case of non-abelian geometric phases. 

Although many of the details have not been worked out for SU{n) groups, 
(the Euler angle parameters, the adjoint representation, etc.) the method of 
identifying the space of the parameters is the same (see 0]). For a system 
with n states, one may express a general diagonal density matrix in terms of 
the squared elements of an n — 1 sphere. Then to take it to a general basis, 
one acts with the appropriate SU(n) matrix. The result is always a subset of 
SU(n)/T n ~ 1 , where T n ~ l is the maximal (n — 1) torus for the group. If there 
are degenerate eigenvalues in the matrix, this space is reduced. For example 
in the case of three states discussed shortly, the parameter space is a subset 
of SU(n)/(SU(2) x U(l)) since there exists a two-fold degeneracy. For an 
m-fold degeneracy we reduce the space by SU{m). In the case of an adiabatic 
approximation, we will see this is a proper subset, but were we to relax this 



condition, the space would be isomorphic to these spaces, not subsets. In 
this way, one may identify a necessary condition for non-abelian geometric 
phases, namely the existence of the degeneracy and thus an SU(m) factor in 
the denominator of the above coset expression. 

Using this parameterization one gains essentially nothing over the ex- 
pression of the Bloch sphere for two-state systems. In that case the common 
parameterization of the Bloch sphere, 

a 
l-o 

is really no different than the one presented here, 

cos 2 9 
sin 2 9 

except that positivity is automatic. However in the case of three state sys- 
tems, we have 

cos 2 # sin 2 

sin 2 # sin 2 ] (9) 

cos 2 . 

This is a convenient parameterization since the analogous Bloch sphere would 
have parameters with a non-rectangular domain. The parameterization given 
here (see also |15|]) then helps with the analysis of three state density matrices 
and their corresponding entropy [pl| . 



4 Connection, Curvature, and Abelian Geo- 
metric Phases 

In the spirit and notation of Nakahara ||, we can now derive the connec- 
tion one form, the curvature and the Geometric Phase of the three state 
system. The connection one form, sometimes called Berry's connection, can 
be written in terms of ip in the following way. Define the total phase to be 



if = i / A 

A = A^dx fl = -i(^\d\?p), (10) 



where d is the ordinary exterior derivative. Using equation (7), this becomes 
A = dx + sin 2 6»[cos 2 (3{da + d>y) - sin 2 j3(da - dj)]. (11) 



This agrees with reference [J 1] if the following identifications are made with 



those quantities on the left being those of reference |Tl| and those on the 
right being ours. 

r] <-> x 6^6 X x <"► a + 7 X 2 ^"-7- 
The corresponding curvature two form is given by 

F = dA= -idift A dip 

= sin 26 cos 2 (3d0 A d(a + 7) - sin 2 6 sin 2(3d(3 A d{a + 7) 

- sin 26 sin 2 (3d6 A d(a - 7) - sin 2 6 sin 2(3d(3 A d{a - 7). (12) 

This then, is the analogue of the "solid angle formula" for the two state 
systems. In other words, the integral of this curvature two form gives the 
geometric phase, just as 

in two state systems, where Q is the solid angle for the two sphere. The 
geometric phase is just the integral of the connection one form without the 
overall phase factor x, that is, 

(f g — sin 2 #[cos 2 (3 {da + dj) — sin 2 (3 {da — d^)} 

= /[sin 2 6 cos2(3da + sin 2 6d-y], (13) 



which again, agrees with [11 



5 Non-abelian Geometric Phases 

In this section a novel way of obtaining geometric phases for 3-state systems 
is given. This method is a generalization and simplification over the method 



presented in QT3 1 and a generalization over the method given in the previous 
section. The way the connection one-forms for the 3-state systems are derived 
here uses the fact that the state space of the system can be expressed in 

8 



terms of the group SU(3). This enables the calculation of the forms without 
diagonalization of the Hamiltonian. In effect, the Hamiltonian is taken to be 
in diagonal form initially. It is then "undiagonalized" by an SU(3) action 
which takes it into a general non-diagonal hermitian matrix. This method 
has the advantage of being potentially generalizable to other states, not just 
eigenstates of the Hamiltonian. (Of course, one has to be careful of what the 
adiabatic assumption means then. This is well described in fT3j .) It also has 
the advantage of being generalizable to SU(n). Whereas one does not have 
a way of finding the eigenvalues of an n x n matrix, one would be able to 
use SU(n) matrices and derive the connection forms for an n-state system. 
(Again, see jI5| .) 

The aim is to find the adiabatic non-abelian geometric phase associated 
to the two-fold degeneracy of energy eigenvalues of the general Hamiltonian 
for a 3-state system. These are the simplest non-abelian geometric phases. 

Let H(t) = H(R(t)) be the time dependent Hamiltonian of the system 
and let E n (t) be its eigenvalues. Then if the Hamiltonian is periodic in time 
with period T, i.e., the curve C:[0, T] — > M is closed. Here M is the manifold 
parameterized by the coordinates R. For the adiabatic approximation, n 
labels the eigenstates, \ip), of the Hamiltonian and does not change. This 
means there is a unitary matrix U(n) relating \ip(T)) and |^(0)) which is 
given by 



-t/c ;En(t) v 



i§ c A n 



Here V is the path-ordering operator and A n is a Lie algebra valued (con- 
nection) one-form whose matrix elements are locally given by: 

A^ = i(n,a,R\d\n,b,R}. (14) 

It is important to note that the Hamiltonian is a 3 x 3 Hermitian matrix 
which can be viewed as an element of the algebra of SU(3), i.e., 

8 

H(R) =b^R i \ i , 

where R % are real parameters, the \ are Ao = I3X3 and the Gell-Mann 
matrices of Table (1). Here the constant b is taken to be one. The adiabaticity 
assumption may then be expressed as T y> 1. 



The Hamiltonian, H, can be expressed in terms of the diagonalized Hamil- 
tonian, Hp. 

H(R) = U(R)H D U- 1 (R), 

where U(R) E SU(3) and 

E x 
H ■ | Ex 
E 3 

In this form it is obvious that M C CP 2 and what is more, it is clear from 
( [|16|1 ) that only the angles a, (3, 7 and 9 will remain since Ai,A2,A3 and X s 
commute with Hjj. Explicitly, the Hamiltonian in undiagonalized form, H, 
is given by 

H n = E ± (cos 2 (3 cos 2 9 + sin 2 /3) + £ 3 cos 2 (3 sin 2 

# 12 = ( J E 1 - J E 3 )e" 2M cos/5sin/3sin 2 ^ 

# 13 = (^a-S^e-^ ^) cos^ sin0 cos 

# 21 = ( J Bi- J B 3 )e 2ia cos / 5sin/5sin 2 e 

# 22 = Ei (sin 2 f3 cos 2 + cos 2 (3) + £ 3 sin 2 /3 sin 2 9 

H 23 = (E 1 -E 3 )e l{a -^sm(3sm9cos9 

H 31 = (E 3 -E 1 )e i{a+ ' y) cos(3sm9cos9 

H 32 = (E 1 -E 3 )e- i{a -^sm(3sm9cos9 

H 33 = E 1 cos 2 9 + E 3 sin 2 9 

It can easily be shown that these angles parameterize CP 2 . In this way 
one can easily identify the patches needed for certain circumstances. This is 
analogous to the calculation here. 

As is well known, the matrix that diagonalizes H is composed of its 
eigenvectors. Therefore, given that H = UHdU^ 1 , H d = U~ X H\J ", so we 
have our \ip)s, the eigenvectors of H, they are 

e -i(a+7) cos p cos g \ / e -*(«-r)sin/3 \ / e -*(<*+7) cos (3 sin 9 

_ e *(a-7) sin p cos 9 , e^ a+7 )cos/3 , -e i( - a ~^ sin(3sm9 

-sin 9 J \ / \ cosfl 

One can check that these are already orthonormal due to the fact that U G 
SU(3). 

10 



Now all that needs to be done is calculate the connection forms given by 
fll4]) . These are given by 

A\ = cos 2(3 cos 2 9 da + cos 2 9 d'y, 

and 

A ( — cos 2(3 da — d r y e~ 2n [sin 2(3 sin 9 da — i sin 9 d(3] 

2 = y e 2 * 7 [sin 2/3 sin 9 da + i sin 9 d(3] cos 2/? sin 2 9 da + sin 2 d 7 

This is a expression in terms of SU(3) Euler angle coordinates. We can 
generalize this by using the expression @. This allows us to express the 
density matrix for an n-state system in terms of the Euler angle coordinates 
and the components of the n — 1 sphere along the diagonal and an overall 
scale factor. Thus the eigenvalues need not be those of the Hamiltonian 
but of any observable. Then a similar analysis holds for states that are not 
eigenvectors of the Hamiltonian but eigenvectors of another observable with 
the caution that, as stated before, one must be careful of what one means by 
an adiabatic approximation. 

6 Conclusions/Comments 

The diagonalized density matrices can be parameterized by the squared el- 
ements of the sphere (an n — 1 sphere for a system of n states) combined 
with an SU(n) action. This novel parameterization helps to identify the pa- 
rameter spaces of these systems. The spaces are isomorphic to subspaces of 
SU{n)/T n for all eigenvalues unique, to subspaces of SU(n)/(SU(2) x T n ~ l ) 
for one two-fold degeneracy, SU{n)/{SU{3) x T n ~ 2 ) for one three-fold de- 
generacy, etc. This is because the density matrix and Hamitonian are both 
in the algebra of the group and can represented as UAU', where U G SU(n) 
and A is the non-diagonal density matrix or Hamiltonian. When this is the 
Hamiltonian, we immediately know the eigenvectors because they are the 
rows of the matrix that diagonalizes the Hamiltonian. This enables the eval- 
uation (in principle) of the geometric phases for the n-state systems. Here 
we have shown this explicitly for the case of three quantum states. 

In this analysis the Euler angle parameterization has been extremely use- 
ful and although its generalization to SU(n) is possible, the decomposition 
into components of the spheres and SU(n) actions is independent of the 
parameterization. 

11 



In |[[B) applications to multi-pole Hamiltonians were discussed. I would 



like to add that there are phenomenological nuclear physics models that 
use SU(3). These multi-pole Hamiltonians are expressible in terms of the 



differential operators in |I6|| , and [17|. The author expects to perform a 



further analysis of the relations to those and other multi-pole Hamiltonians 
in the near future. 
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